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Abstract 



Oh; 

Lattice QCD results reveal that the critical parameters and the order of the quark-hadron 
phase transition are quite sensitive to the number of dynamical flavours and their masses 
included in the theory. Motivated by this result we develop a phenomenological equa- 
tion of state for the quark-gluon plasma consisting of n/ flavours retaining the entropy 



per baryon ratio continuous across the quark-hadron phase boundary. We thus obtain a 
generalised expression for the temperature and baryon chemical potential dependent bag 
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. constant. The results are shown for the realistic case, i.e., involving u, d and s quarks 

^ \ only. We then obtain a phase boundary for an isentropic quark-hadron phase transition 

P-t. 

q h " using Gibbs' criteria. Similarly another phase boundary is obtained for the transition to 

<D : 

an ideal QGP from the solution of the condition -B(/i, T) = 0. The variation of critical 
temperature T c with the number of flavours included in the theory. Also the variation of 
(e — 4P)/T 4 with temperature are studied and compared with lattice results. Finally the 
strange particle ratios 4, § and ^= are obtained at both phase boundaries. We propose 
that their variations with the temperature and baryon chemical potential can be used in 
identifying the quark-gluon plasma in the recent as well as in future heavy-ion experiments. 
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1 Introduction 



Recent lattice gauge simulations reveal interesting results [1-3] for the number of flavours 
and their masses included in the theory, in determining the order and the critical pa- 
rameters of the phase transition. In realistic cases, i.e., for three or more flavours, the 
transition is found to be a first order. However, for two light flavours, one finds a contin- 
uous transition. As the mass of the strange quark is increased from zero, the transition 
changes from first-order to a continuous one. Similarly the critical parameters are also 
sensitive to the number of flavours rif and it is shown that the critical temperature T c 
drops with increasing flavours (T c oc rif~ l l 2 ). The lattice approach, however, is still not a 
suitable tool for studying the strongly interacting matter in the baryon rich environment. 
Therefore, the critical behaviour of such a matter with finite baryon chemical potential 
can reliably be studied in the framework of a QCD-motivated phenomenological models. 
However, within the framework of such models, the specific entropy per baryon (S/B) 
is found to be discontinuous across the phase boundary [4-6]. Leonidov et al. [7] have 
recently proposed a modified equation of state (EOS) for the quark-gluon plasma (QGP) 
with a temperature T and baryon chemical potential \i dependent bag constant in order 
to ensure that the S/ B ratio is continuous along the phase boundary. In our earlier works 
[8-9], we have modified the /i and T dependence of B(/i,T) by incorporating the QCD 
perturbative corrections in the EOS for QGP and we have also explored in detail the 
consequences of such a bag constant on the deconfining phase transition. The purpose 
of this paper is to generalise the expression of B(yU,T) for n/ flavours in the low baryon 
density region and then we investigate the effect of a massless as well as a massive s-quark 
on the critical behaviour of the phase transition. We compare our results with those of 
lattice in the high T region qualitatively. It should be noted that when S/B ratio is 
made continuous across the phase boundary, the quark-hadron phase transition is still 
first-order and it occurs at a common temperature and chemical potential. 
We have then determined the critical parameters from the two phase boundaries - one 
obtained from the Gibbs criteria for an isentropic, equilibrium phase transition and the 



2 



other for the transition to an ideal, noninteracting QGP. Finally we study the variations 
of the ratios ^, § and ^= either with the temperature or with baryon chemical potential 
at both phase boundaries. Here we mainly concentrate on the particle production from 
the QGP in the midrapidity region so that it can be used for studying the properties of 
QGP produced at future RHIC and LHC experiments. 



2 Formulation of B(/i,T) for nj Flavours 

The chemical potential (/x) and temperature dependence (T) of the bag constant can be 
derived by constructing an isentropic, equilibrium phase transition from a quark-gluon 
plasma (QGP) to a hadron gas (HG) at a fixed T and /x, and we put a constraint that 
entropy per baryon ratio is continuous across the phase boundary i.e., 

= | (i) 

QGP -D HG 

The above condition can be achieved by assigning a T and /x dependence to the bag 
constant explicitly [7]. The resulting correction factors — 9B{ ^ T " > and — 9B ^ T ^ to the 
entropy density and the baryon density, respectively, in the QGP phase will modify eq.(l) 
as follows : 
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The above partial differential equation can be solved iteratively by exploiting the symme- 
try property of the bag pressure, i.e., B(/j,,T) = B(— /x, —T). 

We use simple models to describe the EOS in each of the two phases and then perform a 
Maxwell construction to determine their common phase boundary. The QGP phase con- 
sists of a perturbatively interacting gas of quarks and gluons. Hence the thermodynamical 
potential for a QGP can be written in terms of the grand partition function as[10] 

n(T,V,fi) = -^\nZ QGP (3) 



and In Zq G p = \nZ°Q G p + \nZ Int Q G p + \nZ Vac Q G p, where Z° is the zeroth order con- 
tribution, Z Int arises due to perturbative corrections and the term Z Vac represents the 
non-perturbative vacuum contribution in the form of a /i and T dependent bag constant 
-B(/VT). 

The zeroth order contribution to the partition function is : 

^--fe&FjtgL-FfoT^ (4) 
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and the perturbative interaction part up to second order in strong coupling constant g s 
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(a s = ^-) can be written as follows by using finite temperature field theory : 
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(5) 



where Ff =F^ +Fj ; F^(Fj ) are the fermionic (antifermionic) equilibrium distribution 
functions of f-th flavour, respectively : 



F+ = [e Wt-rt + i 



(6) 



e £(%+/*/) + i 



(7) 



Here -E 1 / is the energy of f-flavour. Considering massless quarks the total thermody- 
namic potential for the QGP phase reduces to the following expression [11] : 
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where g(fi f ,T) 

If we further make use of the various conservation conditions such as strangeness (/x s = 
0) and charm conservations, and also put /i u — fx^ — /i q — f, we get the expression for 
the pressure in QGP : 
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(9) 



In these calculations, /x and T dependence of a s is taken as follows [12] 

An 
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where 
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Here /g is the gluonic and Ff is the quark distribution functions, gf is the quark degen- 
eracy factor for the f-th flavour. 

Hadron gas consists of all strange as well as nonstrange hadrons. The partition function 
for baryons (B) then becomes 
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(11) 

here c = 1 - f s and d = 1 - 2f s and / s (= ^) is a factor arising from the strangeness 
conservation condition n s - = in the hadronic phase and n±, ri2 and n 3 are the 
number of nonstrange, singly and doubly strange baryons, respectively. The summation 
arises due to quantum statistics used for the distribution functions. Similarly the partition 
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function for the mesonic (M) sector is 

+ f:^^^(^?)coshA (12) 

n=l 

Here the first term is for massless pions, second and third term is due to nonstrange and 
strange mesons, respectively. In order to incorporate repulsive interactions among the 
hadrons in the HG phase within the excluded-volume approach, we multiply all pointlike 
quantities by a volume correction factor [13] of the type [ 1 + V Q n° HG (T ,/x,)]' 1 where 
n° HG (T,[x) is the net baryon density for pointlike baryons. 

The requirement of conserved entropy per baryon across the phase boundary leads to the 
following expression in the low baryon density and high temperature limit : 

a» - nQGP S m 

We have neglected the terms proportional to fx 2 and fx 3 as well as |^ term. We have 
also assumed that the entropy in the hadron gas at a large T and very small fx is mainly 
contributed by mesons and thus Shg — Sm- 

In the low baryon density limit, we can approximate the entropy density Sq GP as almost 
independent of fx and hence we get, 



- n °QGP ~ n HG (13) 



P B (T,fx)-P B (T,fx~0) 



(14) 



Here P B is the baryonic pressure in the HG phase. Although Sq GP may not be fx- 
independent but it involves fx 2 and higher order terms, which we ignore in the fx — > 
limit. Finally after second iteration, we get the expression for B(/i,T) at high T and small 
fx limit as : 
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Using these equations of state (EOS) for the QGP and HG, one can construct a first-order 
isentropic and equilibrium phase transition with the help of Gibbs criteria i.e., Pqgp = 
Pug and thus determine the critical parameters T c , \i c for the phase boundary. Similarly 
one can also determine the other phase boundary by determining T c /, /v for the phase 
transition from a hadron gas to an ideal QGP from the condition B(/z,T) = in the 
eq.(15). We can then calculate the strange particle ratios along these boundaries. We 
choose T c i as the approximate temperature where (e — 3P)/T 4 = 0, signifying the onset 
of the ideal QGP behaviour (e and P being the energy density and pressure in the QGP 
phase, respectively). Thus we get two sets of critical parameters (T c , /i c ) and (T c /, /v) 
along which we will calculate strange particle ratios. The region under the curve defined 
by (T c , fi c ) is considered effectively as the hadron gas phase, while the ideal QGP phase 
resides above the curve defined by (T c /, /v). The region between the two curves may be 
considered as the transition region from a deconfined phase consisting of massive quarks 
into a chiral symmetry restored ideal QGP phase. The strange particle ratios along (T c , 
fi c ) curve predict quantitavely the minimum QGP values because they represent the de- 
confining phase boundary. Similarly (T c /, /v) curve ensures that, if measurements of the 
observables yield an experimental value for temperature and quark chemical potential 
above this curve, an ideal quark-gluon plasma has been produced in the heavy-ion colli- 
sions. 



3 Calculation of Strange particle Ratios 

The idea that the strangeness abundance can provide a useful signature for quark-gluon 
plasma formation originates from the enhancement of the strange-antistrange quark (ss) 
pairs in a thermally equilibrated and baryon dense QGP as compared to the production of 
u, d light quark flavours[14]. This is possible only when the Fermi energy of the u, d quarks 



become larger than the strange quark mass and hence the Pauli blocking stops further 
creation of light quark pairs. For a qq symmetric and baryonless QGP, the abundance 
of ss can still occur provided the temperature is much larger than the strange quark 
mass. Furthermore, the strangeness abundance in QGP is also related with the lower 
mass threshold involved in the creation of an ss pair in QGP rather than KK pair in the 
hadron gas. Thus it seems worthwhile to investigate in detail the utility of strangeness as 
a QGP signal in the context of future LHC and RHIC experiments where the production 
of baryonless QGP is expected. However, in this context, we will be searching for the 
strange particle ratios ^, § and ^= as a QGP signal. 

In calculating the above ratios we entirely dwell on the formalism in the quark-gluon 
plasma phase formed at a high temperature and a small baryon density. Thus the com- 
position of the QGP mainly involves the gluons and hence the glue-based processes dom- 
inantly contribute to the strangeness production. The probability for the creation of 
quarks of different flavours is calculated in the particle production mechanism [15,16] : 

f q = fo e~— (16) 

where m u = 5 MeV, m d = 9 MeV, m s = 170 MeV, and the QCD string tension k — 1 
(GeV/fm). The normalisation constant f Q is given by 

fo = l — (17) 

The hadronisation of QGP is still not understood properly. We assume that once QGP 
has been produced, the resulting hadrons do not get sufficient time to achieve thermal 
and chemical equilibrium. Thus the chemical composition existing during the QGP stage 
does not change and hence the strangeness abundance achieved in the QGP phase still 
survives during the process of hadronisation. Thus the probability for the creation of a 
primordial hadron with q number of quarks per unit of phase space volume is obtained 

P = H fq\9q l q e"^ (18) 
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where f q is the gluon fragmentation probability, and involves the quark mass dependence, 
the quark fugacity X q = e^l T , g q is the statistical degeneracy factor, j q is the relative equi- 
libration factor and E q is the energy of the q-th quark. We assume complete equilibration 
for u, d, s quarks in QGP, i.e, 7 U = 7d = 7s = 1- The energy of the emitted particle is 
E = J2 q Eq- This factor is integrated out in eq.(18) when we carry out the phase space 
integration and take the particle ratios in the same mr range : E = uit cosh(y pr — y mr ), 



where tjit = \Jm? + °p\ [18], with y pr (y mr ) being the projectile (mid) rapidity. Thus for 
ratios like ^, § and |t=, the numerator as well as denominator gives similar kinematical 
factor after integrations and hence they are cancelled out. However, these ratios depend 
upon the product of f q \ q and hence on the mass of the quark as well as on the tempera- 
ture and quark chemical potential. We introduce A; = e^^ T as the fugacity of i-th hadron 
species in HG and is simply the product of the fugacities of the constituent quarks so that 
A at = X q , \k = X q Xs, etc. Since the u, d, s flavours are separately conserved in the time 
scale of hadronic collisions, the production (or annihilation) can only occur in pairs. So 
the chemical potentials for particle and antiparticle are opposite to each other, we get X q 
= X^ 1 . Finally, particle ratios in the same mx range, take the form [17] : 

A = e -2{^ u+l i d )/T e -2 Ms /T ^ 
| = e -WT e -WT (2Q) 
^_ = e -2^/T e -2 Ms /T ( 21 ) 

All the above ratios depend on the factor exp(±/i q /T) and strange quark chemical poten- 
tial (fi s ) explicitly. In the pure QGP phase, however, /i s is identical to zero because of the 
exact strangeness conservation. The factor 7 S accounts for much of our ignorance about 
the dynamics of strangeness formation and the aproach to equilibration of the strange 
quarks in the QGP phase. A value 7 S = 1 is believed to favour QGP interpretation of the 
data. It is convenient to denote \i q = (/i u + /v)/2; Sfj, = fid — fi u , = 3 \i q ; where fi q is 
quark chemical potential, he is the baryochemical potential and Sfi describes the (small) 
asymmetry in the number of up and down quarks due to neutron excess in heavy-ion 
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collisions, i.e., -| ^ 0.5. In practice we find S/j, to be very small [19], and, therefore, we 
neglect it here. Finally we can calculate these ratios on the phase boundaries using the 
values of the critical parameters (T c , fj, c ) and (T c /, ji c i). 

4 Results and Discussions 

In Fig.l we have shown the variation of B(fi,T) with temperature T as obtained from 
eq.(15) in the low fi and high T limit. We find that the values for B(/z,T) differ much 
from the case of u, d massless quarks and decreases faster as T increases. Moreover, the 
incorporation of interactions into the calculation modifies the result significantly. We find 
that B(/i,T) becomes zero for T = 210 MeV if the interactions among quarks and gluons 
are incorporated and it is far less than the value T = 440 MeV obtained in the case for 
three massless quarks included in the QGP without any interactions as shown by curve 
B. Vanishing of B(yU,T) at a certain value of /z and T signifies that the quarks and gluons 
are almost free. We infer that the ideal gas limit for a QGP consisting of three flavours is 
reached at a value T»T C with the critical temperature T c ~ 160 MeV for the deconfining 
phase transition. However, we find that the calculation shows very little change when we 
take u, d massless quarks and s as a massive quark (m s = 150 MeV) if compared with 
the case of three massless quarks. 

In order to illustrate the above point more clearly, we have plotted in Fig. 2 the variation 
of the quantity (e — 3P)/T 4 with temperature (T) for a QGP consisting of u, d, s mass- 
less quarks. Here e represents the energy density and P is the pressure of the QGP. The 
quantity (e — 3P) represents a measure of the ideal gas behaviour since it vanishes for an 
ideal gas. We find that (e — 3P)/T 4 asymptotically vanishes. This curve agrees with that 
obtained in the lattice gauge results. Thus it gives us confidence in the QCD-motivated 
calculations for the EOS of a QGP. 

We conclude that the inclusion of n/ flavours in the temperature and baryon chemical 
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potential dependent bag constant yields results in agreement with those obtained in the 
lattice gauge results [3]. However, we find that there is very little difference in the results 
whether the strange quark is massless or massive. This result differs from the lattice 
finding. Moreover, our phenomenological model is valid for a QGP with a finite but small 
baryon chemical potential. The B(/x,T) thus obtained guarantees the phase transition be- 
tween QGP and hadron gas at the same value of temperature and chemical potential and 
hence the process of reheating is not required during the mixed phase of hadronisation. 
It also yields continuity of entropy per baryon ratio across the phase boundary. We thus 
hope that the phenomenological form for B(/i,T) thus obtained will be of considerable 
use in deriving the properties and the signals of QGP to be produced at proposed RHIC 
and LHC experiments because we expect the net baryon density in these experiments to 
be low enough as used in the above calculation. However, the baryon density will never 
vanish precisely and, therefore, the lattice calculations can still not be used as such here. 
We, therefore, hope that the description of a deconfinement phase transition achieved in 
a QCD-motivated model by making the bag constant dependent on ji and T will be a 
more realistic one in the future LHC and RHIC experimental situations. 

In Fig. 3, we have separately shown the variations of ^ with the baryon chemical potential 
and also with the temperature. We have separated the regions of HG phase, interacting 
QGP phase and the ideal or non-interacting QGP phase. Our results are valid in the low 
baryon density limit and hence cannot be extended beyond \x > 300 MeV which corre- 
sponds to fi q = 100 MeV. In future colliders like RHIC and LHC, we expect the formation 
of an ideal quark matter at /i < 100 MeV. It means that the value of our calculated ratio 
j > 0.56 should indicate the formation of such a matter. However, the variation of ^ 
with the temperature does not show any similar constraint on the temperature. It simply 
tells us that an ideal QGP can be formed at temperatures T < 180 MeV. The present 
experimental value [20,21] of ^ is 0.22 ± 0.01 and it signifies a hot and dense HG rather 
than a QGP unless \i > 300 MeV for which our predictions are rather unreliable and 
inconclusive. 
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In Fig. 4, we have shown the similar variations for the ratio § with [i and T separately. 
We again find that an extraordinary large value for this ratio ( > 0.74) at pi < 100 MeV 
will indicate the formation of QGP at RHIC or LHC experiments. In Fig. 5, we have 
also shown the variations of ^= either with /i or with T. In this case, we find that the 
formation of an ideal QGP at /i < 100 MeV requires the value of the ratio ^= < 1.3. 

Recently Asprouli and Panagiotou [17] performed an identical analysis. However, they 
fixed the phase boundaries by choosing T c /, /v etc. quite arbitrarily Thus their compar- 
ison with the experimental data revealed a faulty conclusion that the present CERN and 
AGS heavy-ion experimental data [20,21] indicate QGP formation. Furthermore, they 
conclude that Sulphur induced reactions at midrapidity region reveal that the ideal QGP 
has been produced to within 60 % possibility. Our results, on the contrary, suggest that 
the CERN experimental datas correspond to an ideal thermalised HG picture. However, 
the results obtained here clearly indicate that strange particle ratios can provide a signal 
in the midrapidity and baryon-free region. The experimental results do not correspond 
to this region. Nevertheless, the EOS employed here for the QGP phase breaks down 
strictly at n#=0 region because entropy per baryon ratio becomes meaningless. So our 
analysis is applicable only in the low baryon density (ng) region. 

It has been claimed that thermal gluon decay into a quark-antiquark pair dominates for 
a wide range of quark masses in the midrapidity region and thus provides the most im- 
portant source for the quark or the antiquark density [22,23]. Normally, the gluon cannot 
decay into a strange quark-antiquark pair because its thermal mass is below the threshold 
required for such a pair creation. For a temperature (T) around 200 MeV and the coupling 
constant g = 2 for a QGP consisting of two massless flavours, one gets the thermal gluon 
mass m g = | g T = 267 MeV in the lowest order perturbation theory and threshold for ss 
pair creation corresponds to ~ 300 MeV. However, it has been suggested that in addition 
to acquiring a thermal mass of the order of g 2 T, gluons also acquire a width determined 
by the large damping rate [24,25]. Thus g — > ss decay is allowed and thus can account 
for the strangeness enhancement in the midrapidity region. 
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In summary, we have examined the strange particle production in the midrapidity region 
within a QGP formalism. We assume that the gluon fragmentation into qq pairs is the 
main source of the quark and antiquark density in the plasma. However, we have calcu- 
lated the values of the critical parameters for the transition from the HG to an ideal QGP 
using our EOS for the quark matter with a /i and T dependent bag constant obtained 
at a low baryon density. The lower limit to the critical parameters corresponds to Gibbs 
criteria for an equilibrium and isentropic phase transition and upper limit is obtained by 
putting B(fj,,T) = 0. We have plotted the strange particle ratios as a function of T and 
separately along the two phase curves. From the values of these strange particle ratios one 
can infer an approximate values of temperature and the baryon chemical potential reached 
in the heavy-ion experiments. Our studies reveal that the strangeness enhancement can 
still be regarded as the signature for a QGP formed in the central rapidity or at almost 
baryon-free region. However, in the absence of a reliable lattice result for % ^ region, 
our approach is justified. We hope that our results will provide a basis for future studies 
regarding the signals of QGP and their detection in the future collider experiments. 
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Figure Captions 



Fig. 1. Variation of bag constant B(fi,T) with temperature (T) at a baryon chemical 
potential fi = 50 MeV. Curve A represents the free QGP EOS consisting of u, d massless 
flavours. Curve B is for the free u, d, s massless flavours. Curve C stands for u, d massless 
and interacting flavours with QCD scale parameter A = 100 MeV. Curve D represents 
interacting QGP EOS but s-quark is massive ( m s = 150 MeV). Curve E is same as D 
but s-quark is massless. In all these curves, we have used B X J^ = 235 MeV. 

Fig. 2. Variation of (e — 3P)/T 4 which is a measure of an ideal plasma behaviour with 
temperature (T) for /i = 50 MeV and B]^ = 235 MeV. Here e is the energy density and 
P is the pressure of the QGP phase of three massless flavours. 

Fig. 3. Variations of ^ with /i and T are shown separately. The dashed line represents the 
ratio on the phase boundary determined by Gibbs criteria for equilibrium phase transition 
between interacting QGP and HG whereas the solid line denotes the ratio on the phase 
boundary between HG and ideal QGP phases. 

Fig. 4. The notations are the same as Fig. 3 but for the ratio §. 

Fig. 5. The notations are the same as Fig. 3 but for the ratio J^. 
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